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Introduction : 

Hormone  release  from  pituitary  cells  appears  related  to  a 
modification  of  the  basal  electrical  activity  of  the  cell. 
This  electrical  activity  has  been  shown  to  be  modulated  by 
agents  that  also  regulate  secreting  activity.  For  instance 
Thyrotropin-releasing  hormone  (TRH)  triggers  the  release  of 
prolactin  in  GH3  cells  and  simultaneously  leads  to  an 
increase  in  action  potential  frequency  in  these  electrically 
active  cells.  The  membrane  potential  associated  with  the 
appearance  of  released  hormone  in  the  extracellular  fluid 
appears  to  initially  be  hyperpolarized  (thought  to  result  in 

opening  of  a  Ca+  activated  potassium  channel)  followed  by  a 
decrease  in  the  voltage  dependent  K+  currents.  Dubinsky  and 
Oxford  (1)  have  suggested  that  upon  application  of  TRH  1: 
Ca++  is  released  from  intracellular  stores  which  activates 

Ca++  activated  K  channels;  2:  voltage-dependent  K  channel 
openings  are  depressed  during  hyperexcitable  phase  and  that 
3:  TRH  does  not  directly  modulate  calcium  channel  activity. 
During  the  burst  of  action  potentials  during  the 

hyperexcitable  phase,  extracellular  Ca++  enters  the  cells 
through  voltage  gated  Ca++  channels  (perhaps  to  participate 
in  prolactin  secretion)  (2)  and  accumulates  to  a  point  where 
electrical  activity  becomes  again  silent. 

Temperature  effects  the  transition  rates  of  voltage  gated 
channels  (3,4).  This  project  addresses  the  effect  of 
temperature  on  models  of  the  electrical  properties  of 
cellular  action  potentials  and  on  the  accumulation  of 

intracellular  Ca++.  The  work  has  followed  three  directions: 
that  of  developing  computing  tools  to  facilitate  management 
and  display  of  model  results;  developing  a  minimal  model  of 
the  electrical  properties  of  GH3  cells  and  to  develop 
analytic  tools  to  characterize  the  use-and  frequency- 
dependent  properties  of  the  accumulation  intracellular 


calcium  during  bursts  of  action  potentials. 
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Progress  Report: 

1)  Our  research  group  has  been  involved  in  developing 
software  tools  to  facilitate  acquisiton,  management,  analysis 
and  display  of  research  data.  The  most  recent  work  has 
focused  on  developing  tools  for  visualizing  primary  data, 
results  of  simulations  and  displays  of  derived  results. 
Using  the  X-window  system  operating  on  Sun  4  workstations,  we 
have  developed  a  graphical  editor  that  allows  cutting  and 
pasting  segments  of  a  graphic  display.  These  selected 
segments  can  be  printed  or  subjected  to  further  analyses  e.g. 
curve  fitting.  In  addition  we  have  developed  a  tool  for 
scanning  a  sequence  of  research  data  records  (e.g. 
simulations  for  a  set  of  different  conditions) .  These  tools 
have  been  extremely  useful  in  visually  investigating  a  small 
segment  of  a  long  simulation.  The  software  tools  have  been 
described  in  several  manuscripts  listed  at  the  end  of  this 
report  (7,8) . 

2)  GH3  cell  model.  For  a  minimal  model,  we  have  considered  a 
3  component  model:  a  voltage  activated  potassium  current,  1^, 
and  voltage  and  calcium  inactivated  calcium  current,  Ica,  and 
a  leakage  current,  I]_,  that  are  considered  electrically  in 
parallel  with  the  membrane  capacitance.  In  addition,  we  have 
assumed  a  first  order  sequestration  of  internal  calcium. 
Defining  the  membrane  capacitance  as  Cm  we  define  the  minimal 
model  as 


d v 
dt 


dca  +  Ik  +  Il> 

um 


where  V  is  the  membrane  potential,  Ica  is  the  calcium  current 
described  by 


lCa  =  Ica  d(V)  f  <V,Cai)  (V  -  Vca) 

where  d  and  f  are  activation  and  inactivation  gating 
variables,  Vca  is  the  calcium  reversal  potential  and  gca  is 
the  maximum  calcium  conductance.  Similarly,  for  the 
potassium  current 


Ik  -  gkn(V)(V-Vk) 


where  n  (v)  is  the  activation 
potential  and  gk  is  the  maximum 


variable,  Vk  is  the  reversal 
potassium  current.  The  gating 
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variables  are  defined  to  reflect  transitions  in  channel 
protein  conformations  according  to  a  simple  first  order  model 

a 

closed  ^ ^  open 

P 

so  that  for  the  potassium  channel,  n  at  equilibrium  is 
n  =  a/a+p  and  n(t)  is  the  solution  to 

—  =a  (1  -  n)  -p  n 


To  incorporate  Ca++  inactivation  into  the  inactivation 
variable  f,  we  assume  a  first  order  process  leading  to  an 

equilibrium  inactivation  of  the  form  f  oo  [1  +  Ca++/kn]"1 

Finally,  intracellular  calcium  distribution  is  determined  by 


lea  ~  ^ca  ^ai 


where  Ica  represents  calcium  entering  the  cell  through  open 
channels  and  is  normalized  per  unit  volume.  We  have  solved 
these  equations  with  some  preliminary  estimates  of  rate 
constants  and  have  investigate  temperature  dependence 
assuming  a  Qiq  of  3 .  We  have  found  that  this  simple  model 
can  exhibit  3  types  of  temperature  sensitive  behavior  (figure 
1  and  2)  :  simple  oscillations,  bursts  of  oscillations  and 
continuous  oscillations  from  a  depolarized  baseline.  These 
results  suggest  the  nonlinear  terms  in  the  model  can  produce 
behavior  similar  to  that  seen  in  other  systems  exhibiting 
behavior  described  by  the  term,  chaos.  Figure  1  illustrates 
results  when  the  absorption  rate  of  Ca^  is  held  constant 
while  the  temperature  effect  is  restricted  to  channel 
conformation  rates.  Figure  2  illustrates  the  same  sets  of 
rate  constants  but  also  allowing  kca  to  vary  with 
temperature.  The  vertical  axis  is  membrane  potential  (mV) 
while  the  horizontal  axis  is  time  (msec)  .  The  detailed 
mechanism  leading  to  such  dramatic  changes  in  oscillatory 
behavior  is  the  focus  of  current  investigations. 

3)  Our  work  with  analytically  characterizing  the  use- 
dependent  properties  of  Caj_  has  followed  that  of  our  models 
of  ion  channel  blockade  (5).  Basically,  with  each  action 
potential,  intracellular  calcium  is  incremented  by  a 
fraction,  proportional  to  the  difference  between 
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intracellular  and  extracellular  calcium.  For  the  n^h  action 
potential  when  the  channel  is  conducting 


5rci=v(Co  -  ci>  -  kaCl 

where  Cj_  and  CQ  are  intracellular  and  extracellular 

concentrations  and  y  represents  diffusion  rate  down  the 
calcium  concentration  gradient .  When  the  channel  is  not 
conducting,  Cj_  is  reduced  through  intracellular  storage  at  a 
rate  ka  so 


If  the  channel  open  time  is  exponentially  distributed  with 
mean,  tQ,  (5)  then 


Ci  =  C  (oo)  +  [C  (o)  -  C  <oo)  ]  e"  (YC°  +  ka)  to 


where  C(°°)  =  YC0/(y  +  ka) 

While  the  channel  is  not  conducting 

Ci  (t)  =  C  (o)  e"kat 

During  a  burst,  internal  calcium  is  incrementally  increased 
so  that  for  the  n^k  pulse 


•n 


=  C 


ss 


+  (Co  —  Cgg)e 


[kat  r 


+  (y  +  ka)tD]n 


where  tr  is  the  interpulse  interval  between  action  potentials 
during  a  burst  and  tG  is  the  mean  channel  open  time  and 
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p  / 1  0-(Y+ka)t0 

C  (oo)  (1  -e 


-ss 


2  _e-[^r^r  "*■  (Y  ■*"  ka)  1 0] 


Thus,  it  is  possible  to  estimate  the  behavior  of 
intracellular  calcium  during  and  between  bursts  if  the 
channel  open  time  is  exponentially  distributed.  These 
preliminary  analyses  will  be  extended  in  an  attempt  to 
capture  some  of  the  chaotic  features  of  channel  bursting,  and 
to  compare  analytical  results  with  those  derived  from 
numerical  simulation. 


5 


References 


1.  Dubinsky,  J.M.  and  Oxford,  G.S.  Dual  modulation  of 
k  channels  by  thyrotropin-releasing  hormone  in  clonal 
pituitary  cells.  Proc.  Natl.  Acad.  Sci .  82:4282-4286,  1986. 

2.  Benham,  C.D.  Voltage-gated  and  agonist-mediated  rises  in 
intracellular  Ca++  in  rat  clonal  pituitary  cells  (GH3)  held 
under  votage  clamp.  J.  Physiol.  (London)  415:143-158,  1989. 

3.  Hodgkin,  A.L.  and  Huxley,  A.F.  A  quantitative  description 
of  membrane  current  and  its  application  to  conduction  and 
excitation  in  nerve.  J.  Physiol.  (London)  117:500-544, 
1952. 

4.  Partridge,  L.D.  and  Conner,  J.A.  A  mechanism  for  minimizing 
temperature  effects  on  repetitive  fining  frequency.  Amer. 
J.  Physiol.  234 :C155-C161,  1978. 

5.  Starmer,  C.F.  Characterizing  activity  dependent  processes 
with  a  piecewise  exponential  model.  Biometrics,  44:549-559. 
1988. 


6 


Contract  Related  Publications 


1.  Starmer,  C.F.,  Nesterenko/  V.V.,  Undrovinas,  A. I.,  Packer, 
D.L.,  Gilliam,  F.R.,  Grant,  A.O.,  Rosensht raukh,  L.V.  and 
Strauss,  H.C.  Characterizing  ion  channel  blockade  annotated 
files  with  the  guarded  receptor  hypothesis:  Molecular  and 
Cellular  Mechanisms  of  Antiarrhythmic  Agents.  Futura,  Mt . 
Kisco,  NY. 

2.  Starmer,  C.F.  Characterizing  activity  dependent  processes 
with  a  piecewise  exponential  model.  Biometiric  44:549-559, 
1988. 

3.  Starmer,  C.F.,  Undrovinas,  A. I.,  Scamps,  F.,  Vassort,  G., 
Nesterenko,  V.V.  and  Rosenshtraukh,  L.V.  Ethacizin  blockade 

of  Ca++  channels:  A  test  of  the  guarded  receptor 
hypothesis.  Amer.  J.  Physiol,  (in  press). 

4.  Gilliam,  F.R.,  Starmer,  C.F.  and  Grant,  A.O.  Blockade  of 
rabbit  atrial  sodium  channels  by  lidocaine: 
Characterization  of  continuous  and  frequency  dependent 
blocking.  Circ.  Res.  65:723-739,  1989. 

5.  Packer,  D.L.,  Grant,  A.O.  Strauss,  H.C.  and  Starmer,  C.F. 
Characterization  of  concentration  and  use-dependent  effects 
of  quinidine  from  conduction  delay  and  declining  conduction 
velocity  in  canine  Purkinje  fibers.  J.  Clin.  Invest. 
83:2109-2119,  1989. 

6.  Grant,  A.O.,  Dietz,  M.A.,  Gilliam,  F.R.  and  Starmer,  C.F. 
Blockade  of  cardiac  sodium  channels  by  lidocaine  and  its 
hydrophobic  derivative  RAD-242:  Single  channel  analysis. 
Circ.  Res.  (in  press). 

7.  Starmer,  C.F.  and  Dietz,  M.A.  Managing  clinical  research 

data:  Software  tools  for  hypothesis  exploration. 

Environmental  Health  Perspectives  (in  press) . 

8.  Dietz,  M.A.,  Grant,  A.O.  and  Starmer,  C.F.  An  object 
oriented  user  interface  for  analysis  of  biological  data. 
Computers  and  Biomed.  Res.  (in  press). 

9.  Whitcomb,  D.C.  Gilliam,  F.R.,  Starmer,  C.F.  and  Grant,  A.O. 
Marked  QRS  complex  abnormalities  and  sodium  channel 
blockade  by  propoxyphene  reversed  with  lidocaine.  J.  Clin. 
Invest,  (in  press). 


7 


time  (msecs) 

T=17C,  Kca=. 036/msec 


Figure  2 


DISTRIBUTION  LIST 


Stress  Neurochemistry  Program 
Annual,  Final  and  Technical  Reports  (one  copy  each) 

INVESTIGATORS 


Dr.  H.  Elliott  Albers 
Lab.  Neuroendocrin.  &  Behavior 
Depts.  of  Biology  &  Psychology 
Georgia  State  University 
Atlanta,  GA  30303 


Dr.  Gwen  V.  Childs 
Dept,  of  Anatomy  &  Neuroscience 
Univ.  of  Texas  Medical  Branch 
Galveston,  TX  77550 

Dr.  Carl  E.  Creutz 
Dept,  of  Pharmacology 
University  of  Virginia 
Charlottesville,  VA  22908 


Dr.  Mary  F.  Da  11  man 

Dept,  of  Physiology 

University  of  California,  Box  0444 

San  Francisco,  CA  94143-0444 

Dr.  Caleb  E.  Finch 
Dept,  of  Neurobiology 
Univ.  of  Southern  California 
Los  Angeles,  CA  90089-0191 

Dr.  Thackery  S.  Gray 
Department  of  Anatomy 
Loyola  University  Medical  Center 
216  South  First  Avenue 
Maywood,  IL  60153 

Dr.  Richard  F.  Ochillo 
College  of  Pharmacy 
Xavier  Univ.  of  Louisi -ra 
7325  Palmetto  Street 
New  Orleans,  LA  70125 


Dr.  Terry  Reisine 
Dept,  of  Pharmacology 
Univ.  of  Pennsylvania 
School  of  Medicine 
36th  and  Hami 1  ton  Wal k 
Philadelphia,  PA  15104 

Dr.  C.  Frank  Starmer 
P.0.  Box  3181 
Duke  Univ.  Medical  Center 
Durham,  NC  27710 

Dr.  Kent  E.  Vrana 
Dept,  of  Biochemistry 
West  Virginia  School 
of  Medicine 
Morgantown,  WV  26506 


Stress  Neurochemi sty 

Annual,  Final  and  Technical  Reports  (one  copy  each  except  as  noted) 

ADMINISTRATORS 


Scientific  Officer,  Physiology 
Code  1141SB 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217-5000 


Administrator  (2  copies)  (Enclose  OTIC 
Defense  Technical  Information  Center 
Building  5,  Cameron  Station 
Alexandria,  VA  22314 


Program  Manager,  Code  1213 
Human  Factors  Biosciences 
Division 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217-5 000 

Form  50)  Program  Manager,  Code  223 
Support  Technology 
Directorate 

Office  of  Naval  Technology 
800  N.  Quincy  Street 
Arlington,  VA  22217-5000 


Administrative  Contracting  Officer 

ONR  Resident  Representative 

(address  varies  -  obtain  from  business  office) 

Annual  and  Final  Reports  Only  (one  copy  each) 

DoD  ACTIVITIES 


Commanding  Officer 
Naval  Medical  Center 
Washington,  DC  20372 

Commanding  Officer,  Code  404 

Naval  Medical  Research  &  Development  Command 

National  Naval  Medical  Center 

Bethesda,  MD  20814 

Commander 

Chemical  and  8iological  Sciences  Division 
Army  Research  Office,  P.0.  Box  12211 
Research  Triangle  Park,  NC  27709 

Directorate  of  Life  Sciences 

Air  Force  Office  of  Scientific  Research 

Bolling  Air  Force  Base 

Washington,  DC  20332 

Final  and  Technical  Reports  Only 

Director,  Naval  Research  Laboratory  (6  copies) 
Attn:  Technical  Information  Division,  Code  2627 
Washington,  DC  20375 


